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Abstract 

We find the exact winding number distribution of Riemann-Liouville fractional Brow- 
nian motion for large times in two dimensions using the propagator of a free particle. The 
distribution is similar to the Brownian motion case and it is of Cauchy type. In addi- 
tion we find the winding number distribution of fractal time process, i.e. time fractional 



C , Fokker-Planck equation, in the presence of finite size winding center. 
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Fractional diffusion equations are the basic methods to describe a large class of non-equilibrium 
phenomena which show a power law mean square displacement. In the normal case the mean 
square displacement is asymptotically linear in time, i.e. < x^ >~ t, but for the anomalous 
cases it has a non-linear behavior. There are many examples showing anomalous displace- 
ment e.g: protein dynamics [1], relaxation processes and reaction kinetics of proteins |i2j, two 
dimensional rotating flows [3J, porous glasses [1], intercellular transport [5j and so on. There 
are many different kinds of fractional processes but only two of them have more application 
and so have been studied in more detail: fractional Brownian motion (FBM)[6J and the fractal 
time process with fractional Fokker-Planck equation [7j. Both processes are non-Mar kovian 
and have the power law mean square displacement < x^ >^ t^^ , where H \s a real number, 
but they are fundamentally different [8]. 
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The above processes can describe anomalous diffusion of a particle in time or diffusion 
of a macromolecule, in other words it is possible to look at these processes as models for 
dynamics of macromolecules . We have a fractional process which describes the diffusion of 
a polymer in a specific media. One of the most important characteristics of a polymer is the 
winding number distribution which is the simplest quantity describing the entanglement of 
the macromolecule with point-like molecules or molecules with finite size. 

As a simple way to define winding number let's see the process as a two dimensional 
random walker. A two dimensional random walker that starts from the neighborhood of a 
point in the plane tend to follow a path that wraps around that point, the measure of the 
wrapping is given via the winding angle, which is the angle around the reference point swept 
out by walker. Winding angles of paths are of great interest not only from mathematical 
point of view but also because of their application in physics of polymers, flux lines in high 
temperature superconductors and quantum Hall effect [9l[10l[TT]. Winding number of Brow- 
nian motion was calculated long time a go by Spitzer [12]. He showed that the distribution 
of winding number is of Cauchy type, for some other proofs see [131 O HSl [El [17] . Such 
calculation is missing for FBM as a continuous generalization of Brownian motion [6]. In 
this paper we will find the exact asymptotic distribution function of winding number of two 
dimensional FBM in the infinite plane. In addition we will derive the same quantity for the 
fractal time process in the infinite plane with and without finite size winding center. 

The paper is organized as follows: in the next section following the same method as in [18] 
we will find a simple formula similar to Spitzer's result for the winding number distribution 
of FBM. In section 2 we will calculate the same quantity for the fractal time process, in this 
case we will also calculate the winding number distribution in the presence of an obstacle 
which is a good model to describe the entanglement of a polymer in the presence of finite 
size molecule. Finally the last section summarizes our results. 

2 Winding Number Distribution of FBM 

Fractional Brownian motion is defined as a continuous stochastic process with zero mean 
value and < B'jj{t) >~ t^^, with < i/ < 1. To calculate the winding number of such a 
process we will follow the reference [18] which is based on the methods introduced in the 
papers [121 [T3l [T4l [TS] to derive the winding number distribution of Brownian motion. First 
we need to have the probability of finding the particle in r at the time t. This probability is 



calculated for the Riemann-Liouville fractional Brownian motion in [19^ I20j defined as: 

where ^{t) is Gaussian, delta correlated noise < ^{t)S,{t') >= 6{t—t'), and D = frjArni)- This 
process is not markovian and does not have stationary increments, therefore it is different from 
FBM with stationary increments which is more investigated in the mathematical literature. 
The Green function of this process satisfies the following diffusion equation 

^Gj^(ri, r2, t) - ^DH^''-'vl^GH{ri,r2, t) = 6{r2 - ri)6{t), (2.2) 

where V^^ is the Laplace operator acting on ri. The above equation can be considered as the 
effective differential equation describing propagator of FBM (as already proposed in [2 ijjj) . 
Following ^18] one can write the solution of the above equation as 

GHiri,r2,0,t) = -— / exp( ) cosinO) J^{Kr)J^{Kr')KdKdfi, (2.3) 

27r Jo JO ^ 

where J^{s) is the Bessel function of the first kind and ri and r2 are the moduli of ri 
and r2 respectively. It is possible to look at k, as the eigenvalue of the Laplacian then the 
above equation is just a bilinear expansion over the corresponding eigenfunctions. The above 
equation is similar, up to the power of t, to the Brownian motion counterpart and can be 
checked straightly by putting it in to the equation (j2.2p . 

It is worth mentioning exact meaning of the Green function, i.e. G//(ri,r2,^, t), it is the 
statistical weight of trajectories of two dimensional FBM that start at a point ri away from 
the winding center, here origin, and arrive at another point r2 after time t after winding 
angle 9 around origin. If we drop the restriction on the winding around origin then we can 
substitute the integral on /i by the sum and after calculating the sum we will get Gaussian 
distribution function for the trajectories of two dimensional FBM that start at a point ri 
and arrive at point r2 at time t, the variance of this Gaussian function is t , see |19t I20j . 
The reason for including all of the positive /is in the equation ()2.3p is as follows: we are 
trying to calculate winding number distribution in our problem so technically 9 is different 
from 0± 27rn, n is integer, putting sum in the equation (j2.3p will not consider this difference. 
However we should mention that it is also possible to get the true answer by taking the sum 



^it seems that this equation is effectively true and gives the right variance, Green's function and the time- 
dependence of the survival probability [20]. The survival probability is the long time asymptotics of the 
probability Pt that the FBM does not scape from a fixed interval up to time t. 



and following Kholodenco's method for finding winding number distribution [T7]. The other 
point to mention is, if we include also negative /is in the integral then, since J^is) for small 
s is singular for /x < 0, we will not get finite propagator for ri — r2 — > 0. The equation (j2.3p 
is also consistent with propagators in [T71 [191 [2Q| • 

To get winding number distribution we should evaluate the integral (j2.3p for large times. 
The integration over k is truncated at k^ < pi.iu , and so we can replace the Bessel function 
by the first term of its expansion, i.e. J^{s) ~ r(i+u) (§)^- '^^^ integration over k gives 

where Xh = rir2/{2D^t^^). If we choose r2 = ft-'^ + ri then for large t we have \h — JT^P^' 
which goes to zero for large times. If we consider just small Xh in the integrand then the 
integral will be dominated by small /j,. The first order approximation gives 

G„(n,..,^.t)^ J,,„(,^(g'^^, . (2.5) 

To get the distribution for 6, since Xh is small we can use In A// ~ — Ini , after renor- 
malization we have a Cauchy type distribution 

G{x = ^-?— ) = --^, (2.6) 

^ Hint' 7rl + x2' ^ ^ 

for t goes to infinity. For H = ^ this is the same as Spitzer's result. The same result is 



accessible for the other generalization of Brownian motion defined by Bnit) = D J^ [, 
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see [ini, with the same constant D. This process has the same propagator but has a different 
autocorrelation |19t I22j . Since we just need the propagator to find winding number distri- 
bution, the winding distribution of this process is the same as Riemann-Liouville fractional 
Brownian motion. 

It seems that the same calculation should be tractable for different boundary conditions 
if the propagator of FBM could be calculated, of course translational invariance plays crucial 
role in the above argument. The boundary condition corresponding to the problem of walker 
in the presence of finite size winding center is equivalent to solving the equation (j2.2p for 
cylindrical boundary condition with zero generating function on the boundary. We have 
unfortunately not been able to provide a solution for this problem. However it is possible to 
do the calculation for the fractal time process, fractional Fokker-Planck equation, which is 
the main subject of the next section. 



3 Winding Number Distribution of Fractional Time Process 

To define fractional time process or fractal time random walk we will follow the approach 
of [23]. Consider a continuous time random walk so that the waiting time between two 
jumps and the length of the jumps come from a special pdf, ij){x,t). Then one can get the 
jump length pdf, A(x), and waiting time pdf, w(t), just by integrating of il){x,t) on t and x 
respectively. Using 'ip{x,t), the fractal time random walk can be described by the following 
master equation 

/oo roo 

dx' dt'r]{x',t')7p{x-x',t-t')+5{x)5{t), (3.1) 
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where in the process we have taken into account w{t) ~ t^'-^^^^), and the jump's length 
variance is also finite. We are not going to describe all of the interesting properties of 
this process, the only thing we would like to mention is that this process has the following 
fractional Fokker-Planck equation 



-GH{ri,r2,t) - D^-^^^Vl^GH{ruV2, t) = 5{r2 - ri)5{t). (3.2) 

In the above formula we used the Riemann-Liouville fractional derivative ^.i-jj) = r(2H-i') 
/o (f_JYi-iu , with 1/2 < H < 1, see [21]. The corresponding derivative for < if < 1/2 can 
also be defined by adding ordinary derivative to the fractional one. Using the above equation 
it is easy to find the winding number distribution of this process in two dimensions in the 
presence of a disc-like obstacle in the origin which imposing cylindrical boundary condition 
for the above propagator. Following the same steps as in [25J we can find the solution for 
cylindrical boundary condition with zero generating function on the boundary. In the Laplace 
transform space of t, i.e p, we have 

Vl^GHiri,r2,p) - KGH{ri,r2,p) = -a6(r2 -ri), (3.3) 

where k ~ p"^^ , and a ~ p^^~^. If we now go to the Fourier space of 9, i.e fi, we will have 

d \ d o Li ol 

(^-^ H ■^—)GH{ri,r2,p,fi)-{K + ^)GH{ri,r2,p, fi) = 6{ri-r2). (3.4) 

ori ri ori rf ri 

The solution for the above equation is well known, see for example [25j , and has the following 
form 



where we take r2 > ri and R is the radius of disc removed from the plan and the functions 
K^ and I^ are the modified Bessel functions. To get winding number distribution we need 
to integrate over the radial coordinate of the final points. Returning to the 9 and t space by 
inverse Fourier and inverse Laplace transform we can write 

P(M)^/_>.e-^C(. /'''''"'--'^'-'-'";^};^;'"«>^---'''"'-'' ). ,3,6) 

where C{^) = /^°° dr2r2K^{t~^r2), for ri close to the border of obstacle. To get the inverse 
Laplace transform we used steepest descent approximation similar to [25] for large times 
which is like substituting k with t~^ after integrating on p spacqj. The above equation is 
quite the same as the equation (2.15) in [25] with just modified k, so by following Rudnick 
and Hu [25] one can argue that the integrand vanishes exponentially at ±oo for large times. 
Then it is possible to calculate the integral by transforming the integration over /u in to an 
integral around a closed contour in the complex plane as the standard contour integration. 
The integral were calculated in |25j for generic small k, i.e. in the large time limit, and has 
the following form 

^^"^ = m^^ = 4cosh2(2^)- ^^-^^ 

We can not get the limit of zero winding center from the above equation, to get this limit 
we should go back to the original equation of P{0^ t) and take the limit -R ^ 0. It is not 
difficult to see that just the function 1^ in the integral will survive and we will get the similar 
equation as we found for Riemann- Liouville fractional Brownian motion. It is possible also 
to get the same result as ()3.7p by following [18j, which is the same method but with slightly 
different technicality. 

4 Discussion and Summary 

Using the Cauchy distribution that we found for both FBM and fractional time process it 
is not difficult to show that < e >~ tstt. It is the same as the result for FBM with 
autocorrelation < Bnit + T)BH{t) >= ^{{t + Tf^ + t^^ - T^^), which has stationary 
increments, at large n (26j. We should emphasize that our result is not necessarily true for 
this kind of FBM which has different kind of Fokker-Planck equation. However it seems 



For more detail about the inverse Laplace transform see [18] and for the application of steepest descend 
method for inverse Laplace transform see [27] 



that the power law behavior of generating function of winding angle at large times for large 
winding numbers is the general property of different kinds of fractional Brownian processes. 
Finally we think that it is interesting to calculate the same winding distributions for other 
familiar situations specially those proposed in [18j, such as: paths with fixed endpoints with 
and without obstacle in winding center and paths with glued endpoints. These cases as 
argued in [18] are more related to the polymer applications. 

It is interesting also to calculate the winding number distribution by using other methods, 
specially path integral method, and see the connection to fractional quantum mechanics as 
already done for the ordinary Brownian motion in [131 [15] . It is also interesting to verify our 
formula by numerical studies. 

In conclusion we calculated the exact winding number distribution of Riemann-Liouville 
fractional Brownian motion for the point-like winding center which is the first result in this 
context. We did more general calculation for the fractal time process with finite size winding 
center and showed that the equation is, up to a parameter H, similar to the Brownian 
motion winding number distribution. Moreover we showed that at least for large times and 
large winding numbers the generating function of winding number is power law in the absence 
of obstacle. 
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